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ABSTRACT 

Analysis of ten years of high-precision timing data on the millisecond pulsar 
PSR J0437— 4715 has resulted in a model-independent kinematic distance based 
on an apparent orbital period derivative, Pb, determined at the 1.5% level of 
precision (D^ = 157.0±2.4pc), making it one of the most accurate stellar distance 
estimates published to date. The discrepancy between this measurement and a 
previously published parallax distance estimate is attributed to errors in the 
DE200 Solar System ephemerides. The precise measurement of Pb allows a limit 
on the variation of Newton's gravitational constant, \G/G\ < 23 x 10~ 12 yr _1 . We 
also constrain any anomalous acceleration along the line of sight to the pulsar 
to |a /c| < 1.5 x lCT 18 s -1 at 95% confidence, and derive a pulsar mass, m psr = 
1.76 ± 0.20 Mq, one of the highest estimates so far obtained. 

Subject headings: stars: distances — pulsars: individual (PSR J0437— 4715) - 
stars: neutron 



1. Introduction 



In Il993l . IJohnston et al. I reported the discovery of PSR J0437— 4715, the nearest and 
brightest millisecond pulsar known. Within a year, the white dwarf companion and pulsar 
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wind bow shock were obser ved (IBell. Bailes fe Bessell Ill993l ) and pulsed X-rays were detected 
( IBecker fe Triimper Ill993l ). The proper motion and an initial estimate of the parallax were 
later presented alo ng with evidence for se cular change in the inclination angle of the orbit due 



to proper motion (ISandhu et al. 1119971 ). Using high time resolution instrumentation, the 



three-dimensional orbital geom etry of the binary system was determined, enabling a new 
test of general relativity (GR; Ivan Straten et al. 1120011 ). Most recently, mul ti- frequency 



obser vations were used to compute the dispersion measure structure function (I You et al. 



20071 ). quantifying the turbulent character of the interstellar medium towards this pulsar. 



The high proper motion and proximity of PSR J0437— 4715 led to the prediction (IBell fe Bailes 



19961 ) that a distance measurement independent of parallax would be available within a 
decade, when the orbital period derivative (Pb) would be determined to high accuracy. Even 
if the predicted precision of ~1% would not be achieved, such a measurement would be sig- 
nificant given the strong dependence of most methods of distance determination on relatively 
poorly constrained models and the typically large errors on parallax measurements. Even 
for nearby stars, bot h the Hubble Space T elescope and the Hipparcos satellite give typical 



distance errors of 3% (1 Vails- Gabaud 1120071) and so far only t wo distances beyond 100 pc have 
been determined at ~ 1% uncertainty ([Torres et al. 1120071 ). This kinematic distance is one 



of the few model-independent methods that does not rely upon the motion of the Earth 
around the Sun. 



As demonstrated by iDamour fe Taylor I (119911 ). P\> can also be used to constrain the 
variat ion of Newto n's gravitational constant. The best such limit from pulsar timing to 
date JTavlor I Ejjajj , \G/G\ 



(4 ± 5) x 10 



-12 yr -l 



from PSR B1913+16) is c ompromised 



due t o the poorly constrained equation of state for the neutron star companion (INordtvedt 



1990h. The slightly weaker but more reliable limit of \G/G\ = (-9 ± 18) x lO^yr" 1 
(jrCaspi. Taylor fe Ryba 1994 . from PSR B1855+09, which has a white dwarf companion) 
should therefor e be c onsidered instead. A more stringent limit can be obtained from the 

(2.5 ± 3.3) x HT 12 yr- 1 (at 95% 



Splaver et al. I (120051 ) timing of PSR J1713+0747: \G/G\ 
certainty). This limit is, however, based upon the formal errors of Pb, Pb and parallax, which 
are easily underestima ted by standard metho dologies, as we shall demonstrate later. Because 
of this we believe the ISplaver et al. I (120051 ) limit is probably underestimated, but still of 
relevance. However, none of these limits are as strong as that put by lunar laser ranging 



(LLR: IWilliams. Turvshev fe Boggs 1120041 ): G/G = (4 ± 9) x 10- 13 yr -1 . Besides limiting 



alternative theories of gravity, bounds on G can also be used to co nstrain variations of the 



Astro nomical Unit (AU). Current planetary radar experiments (IKrasinsky fe Brumberg 
20041 ) have measured a significant linear increase of dAU/dt = 0.15 ± 0.04m yr -1 , which 
may imply G/G = (—10 ± 3) x 10~ 13 yr -1 , just beyond the sensitivity of the limits listed 
above. 
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As mentioned before, the equation of state for dense neutron star matter is very poorly 
constrained. Pulsar mass determinations can probe the range of permissi ble pulsar masses 
and thereby limit possible equations of state (ILattimer fc Prakash 1 120071 ). Presently, only 
the pulsars NGC 6440 B, Terzan 5 1 and T e rzan 5 J have predicte d masses higher than the 



typical value of 1.4 M^ tjFreire et al. 1 120071 : iRansom et al. 1120051 ); however, as discussed in 



more detail in Section [51 such predictions do not represent objective mass estimates. 

The remainder of this paper is structured as follows: Section [2] describes the observa- 
tions, data analysis and general timing solution for PSR J0437— 4715. Section [3] describes 
how the measurement of Pb leads to a new and highly precise distance. In Section HI this 
measurement is combined with the parallax distance to derive limits on G and the Solar 
System acceleration. Section \5\ presents the newly revised pulsar mass and our conclusions 
are summarised in Section [HI 



Observations and Data Reduction 



Observations of PSR J0437— 4715 were made over a time span of ten years (see Figure 
[1]), using the Parkes 64-m radio tel escope, two 20 cm receiving s ystems (the central beam 



of the Parkes multi-beam receiver (jStaveley-Smith et al. I Il996l ) and the H-OH receiver) 



and four generations of digital instrumentation (see Tabled]): the Fast Pulsar Timing Ma- 
chine (FPTM), the S2 VLBI recorder, and the Caltech-Parkes-Swinburne Recorders (CPSR 
and CPSR2). The FPTM is an autocorrelation spectrometer, whereas the three other in- 
struments are baseband data recording and processing systems that employ phase-coherent 
dispersion removal. 



2.1. Arrival Time Estimation 



For the FPTM, S2, and CPSR backends, the uncalibrated p olarization data were com- 
bined to form the polarimetric invariant interval (IBritton 1120001 ) and each observation was 
integrated in time and frequency before pulse arrival times were calculated through stan- 
dard cross-correlation with an instrument - depen dent template profile. For the CPSR2 data, 
the technique described by Ivan Straten I (120041 ) was used to calibrate 5 days of intensive 
PSR J0437-4715 observations made on 2003 July 19 to 21, 2003 August 29, and 2005 July 
24. The calibrated data were integrated to form a polarimetric template profile with an inte- 
gration length of approximately 40 hours and f requency resolution of 500 kHz. This template 
profile and Matrix Template Matching (MTM, Ivan Straten 1120061 ) were used to calibrate the 
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three years of CPSR2 data. An independent MTM fit was performed on each five-minute 
in tegration, producin g a unique solution in each frequency channel, as shown in Figure 2 
of Ivan Straten I (120061 ) . The calibrated data were then integrated in frequency to produce a 
single full-polarization profile at each epoch. MTM was then used to derive time-of-arrival 
(TOA) estimates from each calibrated, five-minute integration. The application of MTM 
during the calibration and timing stages reduced the weighted rms of the CPSR2 post-fit 
timing residuals by a factor of two. A ll the data reduction described above wa s performed 
using the psrchive software package (IHotan. van Straten fc Manchester 1 120041 ) . 



2.2. Timing Analysis 



Most data were recorded at a wavelength of 20 cm; however, in the final three years, 
simultaneous observations at 10 and 50 cm were used to measure temporal variations of 
the interstellar di spersion delay (co rrections for these variations were implemented in a way 
similar to that of lYou et al. 1120071 ). A linear trend of these delays was also obtained for the 
year of FPTM data, using data at slightly different frequencies close to 1400 MHz. 

The arrival times were analyzed using the Tempo2 pulsar timing soft ware package 
(IHobbs. Edwards fc Manchester 1 120061 : Edwards. Hobbs fc Manchester 1 120061 ) and consis- 
tency with the earlier program, Tempo, was verified. The timing model (see Ta ble [2]) 

and 
The 



is based on the relativistic binary model first derived by iDamour fc Deruelle 



expanded to contain the geometric orbital terms described by iKopeikin I (11995 



(1986) 



1996 



model is optimised through a standard weighted least-squares fit in which all parameters 
are allowed to vary, including the unknown time delays between data from different instru- 
ments, but excluding the mean value of dispersion measure, which is determined from the 
simultaneous CPSR2, 64MHz-wide bands centred at 1341 and 1405 MHz. 

A major differen ce between our implementation of solutions for the orb ital angles Q and i 
and previous efforts (Ivan Straten et al. Il200ll ; IHotan. Bailes &: Ord 1120061 ) is that they were 
implemented as part of the standard fitting routine. This ensures any covariances between 
these and other parameters (most importantly the periastron advance and companion mass, 
see Table [2] and Section [5]) are properly accounted for, thereby yielding a more reliable 
measurement error. The previous works mentioned above derived these effects from an 
independent mapping of % 2 space, leaving the errors of other parameters unaffected. 

As can be seen from Figure (U there are significant low-frequency structures present in 
the timing residual data. Since the standard least-squares fitting routine used in Tempo2 
does not account for the effect of such correlations on parameter estimation, we performed a 
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Monte-Carlo simulation where data sets with a post-fit power spectrum statistically consis- 
tent with that of the PSR J0437— 4715 data were used to determine the parameter estima- 
tions uncertainties in the presence of realistic low frequency noise. These errors, as well as 
the factors by which the original errors were underestimated, are shown in Table [2J As an 
example, the distribution of derived pulsar masses from the Monte-Carlo simulation is given 
in Figure [2j Because of the dispersion measure corrections implemented in the final three 
years of data, one can expect the spectrum of these most precise data points to contain less 
low-frequency noise than the ten year data set as a whole. We therefore expect the errors 
resulting from this analysis to be slightly overestimated. Ongoing research into extending 
the fitting routine with reliable whitening schemes to avoid the spectral leakage and hence 
improve the reliability of the measured parameters, is expected to reduce these errors by fac- 
tors of around two. All errors given in this paper are those resulting from the Monte-Carlo 
simulations, unless otherwise stated. The simulations also showed that any biases resulting 
from the red noise are statistically negligible for the reported parameters. (A full description 
of this Monte-Carlo technique and the whitening schemes mentioned will be detailed in a 
future publication.) 



2.3. Solar System Ephemerides 



Pulsar timing results are dependent on accurate ephemerides for the Solar Sy stem bod- 



ies. The results presented in this paper were obtained using the DE405 model (IStandish 



20041 ) and, for comparison, selected parameters obtained with the earlier DE200 model 
are shown in Table [3j The greatly reduced \ 2 indicates that the newer Solar System 
ephemerides are superior to the earlier DE200, rein forcing; similar conclusions of other authors 
( jSplaver et al. 1120051 ; iHotan. Bailes fc Ord II2006I ). We notice the parallax value changes by 
more than 10 a, and that the different derived values are closely correlated with the ephemeris 
used. Although the effect is not as dramatic as it appears because of the under-estimation 
of the Tempo2 errors, the fact that the DE405 results agree much better with the more 
accurate kinematic distance (discussed in the next Section), strongly suggests that the dif- 
ferences are due to the ephemeris used and confirms that the DE405 ephemeris is superior. 
Finally, we note that the DE405 measurement of Cj (0.016 ± 0.008 °yr _1 ) is consistent with 
the GR prediction for this system (0.0172 ± 0.0009 °yr -1 )- 
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Kinematic Distance 



As shown in Figure [3l the long-term timing history enables precise measurement of 
the orbital period derivative, Pb= (3.73 ± 0.06) x 10~ 12 . This observed value represents a 
combination of phenomena that are intrinsic to the binary system and dynamical effects that 
result in both real a nd apparent accelerations of the binary system along the line of sight 
feell fc Baileslll996h : i.e. 



>obs 



pint 



+ P h Gal + P 



kin 
b 



X 



1 b — 1 h ' b 

where "obs" and "int" refer to the observed and intrinsic values; "Gal" and "kin" are the 
Galactic and kinematic contributions. 

Intrinsic orbital decay is a result of energy loss typically due to effects such as atmo- 
spheric drag and tidal dissipation; however, in a neutron star-white dwarf binary system 
like PSR J0437— 4715, energy loss is dominated by quad rupolar gravitational wave emission. 
For this system, GR predicts (ITavlor & Weisberg Ill982l ) P GR = —4.2 x 10~ 16 , two orders of 



magnitude smaller than the uncertainty in the measured value of P^. 

Galactic contributions to the obs erved orbital period deri vative include differential ro- 
tation and gravitational acceleration (IDamour fc Taylor Ill99ll ). The differential rotation in 
the plane of the Galaxy is estimated from the Galactic longitude of the pulsar and the Galac- 
tocentric distance and circular velocity of the Sun. Acceleration i n the Galactic gravitationa l 
potential varies as a function of height above the Galactic plane (IHolmberg fc Flynn 1120041 ). 
which may be estimated using the parallax distance and the Galactic latitude of the pulsar. 
Combining these terms gives P£ al = (—1.8 — 0.5) x 10~ 14 = —2.3 x 10~ 14 , which is of the 
same order as the current measurement error. 



Given the negligible intrinsic contribution, Equation [T] can be simplified an d rewritten 



in ter ms of the dominant kinematic contribution known as the Shklovskii effect (jShklovskii 



19701 ). an apparent acceleration resulting from the non-linear increase in radial distance as 



the pulsar moves across the plane perpendicular to the line of sight; quantified by the proper 
motion, /i, and distance D from the Earth: 



pobs 



pGal 



P 



kin 



(2) 



where c is the vacuum speed of light. Using t he measured values o f /x, P^, and P^, Equation 
[2] is used to derive the kinematic distance (IBell fc Bailes 1 119961 ): D k = 157.0 ± 2.4 pc. 
This distance is consistent with the one derived from parallax (D n = 150 ±12 pc - see 
also Figure H]) and is, with a relati ye error of 1.5%, com parable in precision to the best 



parallax measurements from VLBI (ITorres et al. 1 120071 ) and better than typical relative 



errors provided by the Hipparcos and Hubble space telescopes (IValls-Gabaud 



2007J) 
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Given the dependence of parallax distances on ephemerides, as described in Section 12.31 
it is interesting to note the robustness of D^. Also, Table [2] shows that the presence of red 
noise corrupts the parallax error by a factor of 7.9, whereas Pb is only affected by a factor 
of 2.5. These facts clearly indicate the higher reliability of as compared to D v . 



4. Limits on P b Anomalies: G and the Acceleration of the Solar System 



Any anomalous orbital period derivative can be constrained by substituting the parallax 
distance into Equation 2, yielding 



p, \ excess 



oGal 



i>kin\ 



(3.2 ±5.7) x l(T 19 s- ] 



(3) 



i n whi ch the error is almost exclusively due to the parallax uncertainty. Following iDamour fc Taylor 
( 119911 ). this can be translated into a limit on the time derivative of Newton's gravitational 
constant (given are 95% confidence levels): 



G 
G 



1 /Ph 



2 VP h 



-5 ±18) x l(T 12 yr -1 



(4) 



This limit is of the same order as those previously derived from pulsar timing (see Section 
[1]), but a currently ongoing VLBI campaign on this pulsar is expected to improve significantly 
on our parallax measu rement, and this should improve our limit to close to that put by LLR 
((4 ±9) x 10~ 13 yr -1 ; IWilliams. Turyshev fc Boggs 112004 ). The LLR experiment is based 



on a complex n-body relativistic model of the planets that incorporates over 140 estimated 
parameters, such as elastic deformation, rotational dissipation and two tidal dissipation 
parameters. In contrast, the PSR J0437— 4715 timing result is dependent on a different set 
of models and assumptions, and therefore provides a useful independent confirmation of the 
LLR result. 

A recent i nvestigation into the possible c auses of a measured variability of the Astronom- 
ical Unit (AU ; iKrasinsky fc Brumberg 1120041) has refuted all but two sources of the measured 
value of dAU jdt = 0.15±0.04m/yr. IKrasinsky fc Brumberg I (120041 ) state that the measured 
linear increase in the AU would be due to either systematic effects or to a time- variation of 
G at the level of G/G = (—10 ± 3) x 10 _13 yr _1 , comparable to, but inconsistent with, the 
LLR limit. 

The anomalous Pb measurements of a number of millisecond pulsars have also been used 
to place limits on the acceleration of the Solar System due to any nearby stars or undetected 
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massive planets (jZakamska fc Tremaine II2005I ). The PSR J0437— 4715 data set limits any 
anomalous Solar System acceleration to |a®/c| < 1.5 x 10~ 18 s _1 in the direction of the 
pulsar with 95% certainty. This rules out any Jupiter-mass planets at distances less than 
117 AU along the line of sight, corresponding to orbital periods of up to 1270 years. Similarly, 
this analysis excludes any Jupi ter-mass planets orbiting PSR J0437— 4715 between ~5 and 
117 AU along the line of sight. IZakamska fc Tremaine I (120051 ) also compared the sensitivity 
of this limit to that of optical and infra-red searches for trans-Neptunian objects (TNOs) 
and concluded that beyond ~ 300 AU the acceleration limit becomes more sensitive than the 
alternative searches. At a distance of 300 AU from the Sun, the 95% confidence upper limit 
on the mass of a possible TNO (in the direction of the pulsar) is 6.8 Jupiter masses. The 
precise VLBI measurement of parallax mentioned above might decrease this to close to one 
Jupiter mass. 



5. Pulsar Mass 

A combination of the mass function and a measurement of the Sha piro delay range can 



be use d to obtain a measurement of the pulsar mass. Using this metho d. Ivan Straten et al 



(2001) derived a mass for PSR J0437-4715 of 1.58 ± 0.18 M® whereas Hot an. Bailes fc Ord 



( 120061 ) obtained 1.3 ± 0.2 M®. It should be noted, however, that these values resulted fr om a 



model that incorporated geometric parameters first described by lKopeikin I (119951 . Il996l ). but 
covariances between these and other timing parameters (most importantly the companion 
mass or Shapiro delay range) were not taken into account. Whilst the length of the data 
sets used by these authors were only a few years, it can also be expected that some spectral 
leakage from low-frequency noise was unaccounted in the errors of these previously published 
values. As described in Section [2J the Monte-Carlo simulations and extended fitting rou- 
tines implemented for the results reported in this paper do include these covariances and 
spectral leakage; it can therefore be claimed that the current estimates (at 68% confidence) 
of m c = 0.254 ± 0.018 M & and m psT = 1.76 ± 0.20 M®, for the white dwarf companion and 
pulsar respectively, reflect the measurement uncertainty more realistically than any previous 
estimate. The distribution of m psr that follows from the 5000 Monte-Carlo realizations is 
shown in Figure [21 together with a Gaussian with mean 1.76 and standard deviation 0.20. 
This demonstrates the symmetric distribution of the pulsar mass likelihood distribution, 
induced by the precise determination of the orbital inclination angle. 

We also note that the new mass measurement of PSR J0437— 4715 is the highest ob- 
tained for any pulsar to date. Distinction needs to be made between the obj ective mass esti- 



mate presented in this paper and the subjective mass predictions presented in lRansom et al. 
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( 120051 ) and iFreire et al. I (120071 ) . The pulsar mass confidence interval presented in this paper 
is derived from the measurement uncertainties of all relevant model parameters, including 
the well-d etermined orbital incl ination angle, i. In contr ast, % is unknown in the Terzan 



51 and J (IRansom et al. 1 120051 ) and PSR J1748-2021B (IFreire et al. 1 120071 ) binary sys- 



tems, and the posterior probability intervals for the pulsar masses presented in these works 
are based upon the prior assumption of a uniform distribution of cos 2. These fundamen- 
tal differences must be accounted for in any subsequent hypothesis testing. Consequently, 
PSR J0437— 4715 is currently the only pulsar to provide reliable const raints on equations of 



state b ased on hyperons and Bose-Einstein condensates as described by lLattimer fc Prakash 



(120071 ). Simulations with Tempo2 indicate that a forthcoming observational campaign with 
a new generation of backend systems can be expected to increase the significance of this 
measurement by another factor of about two in the next year. 



6. Conclusions 

We have presented results from the highest-precision long-term timing campaign to date. 
With an overall residual rms of 199 ns, the 10 years of timing data on PSR J0437— 4715 have 
provided a precise measurement of the orbital period derivative, P^, leading to the first 
accurate kinematic distance to a millisecond pulsar: = 157.0 ± 2.4 pc. Application of 
this method to other pulsars in the future can be expected to improve distance estimates to 



other binary pulsar systems (IBell fc Bailes II1996I ). 



Another analysis based on the P^ measurement places a limit on the temporal variation 
of Newton's gravitational constant. We find a bound comparable to the best so far derived 
from pulsar timing: G/G = (—5 ± 18) x 10 -12 yr -1 . An ongoing VLBI campaign on this 
pulsar is expected to improve this limit, enabling an independent confirmation of the LLR 
limit. 

Previous estimates of the mass of PSR J0437— 4715 have been revised upwards to 
m psr = 1.76 ± 0.20 Mq, which now makes it one of the few pulsars with such a heavy mass 
measurement. A new generation of backend instruments, dedicated observing campaigns, 
and data prewhitening techniques currently under development should decrease the error in 
this measurement enough to significantly rule out various equations of state for dense nuclear 
matter. 
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rms = 0.1 99 {is 




1 996 199S 2000 2002 2004 2006 



Fig. 1. — Combined 20 cm post-fit timing residuals for new and archival PSR J0437— 4715 
timing data. Vertical dashed lines separate the different instruments. 




Fig. 2. — Pulsar mass probability distribution. The solid line shows the histogram of 5000 
pulsar masses derived from a Monte-Carlo simulation with power spectrum and sampling 
equal to that of the PSR J0437— 4715 data set. The dashed line is a Gaussian distribution 
with a mean value of m psr = 1.76 M Q and standard deviation of 0.20 M & . 
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Fig. 3. — Variations in epoch of periastron passage (To) due to apparent orbital period 
increase. A steady increase in orbital period is equivalent to a quadratic increase in T 
relative to periastron times for a constant orbital period. For this plot, T was measured 
on data spans of up to 120 days with a model having no orbital period derivative. The 
formal one-cr measurement errors reported by TEMP02 are shown by vertical error bars and 
the epochs over which the measurements were made are shown by horizontal bars. As the 
mean measurement time was determined through a weighted average of the data contained 
in the fit, these horizontal bars need not be centred at the mid time associated with the 
measurement. The parabola shows the effect of the Pb value obtained from a fit to the data 
shown in Figure [H 
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Day of year 



Fig. 4. — Parallax signature of PSR J0437— 4715. Top: Timing residuals for 
PSR J0437— 4715 as a function of day of year (starting on 18 November), without paral- 
lax but with all remaining parameters at their best-fit values. The smooth curve represents 
the model fit of a parallax of 6.65 mas. Bottom: The same timing residuals with parallax 
included in the model. The overall rms for the top and bottom plots is 524 and 199 ns 
respectively. The double-humped signature specific to parallax originates from the delay in 
pulse time-of-arrival (TOA) as the Earth orbits the Sun and samples different parts of the 
curved wave-front originating at the pulsar. 
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Table 1. Characteristics of the timing data for the four instruments used 



Backend 


Date 


range 


Ref. 


Bandwidth 


RMS 


Observation 


Number of 


TOA 












Residual 


length a 


TOAs 


error a 


FPTM 


1996 Apr - 


- 1997 May 


1,2 


256 MHz 


368 ns 


10 min 


207 


500 ns 


S2 


1997 Jul - 


- 1998 Apr 


3 


16 MHz 


210 ns 


120 min 


117 


160 ns 


CPSR 


1998 Aug 


- 2002 Aug 


3 


20 MHz 


218 ns 


15 min 


1782 


250 ns 


CPSR2 


2002 Nov 


- 2006 Mar 


4 


2 x 64 MHz b 


164 ns 


60 min 


741 


150 ns 



a Displayed are typical values only. 

b CPSR2 records two adjacent 64 MHz bands simultaneously at 20 cm. 



References. (1) ISandhu et al. I (119970 : (2) ISandhu I (1200 lh : (3) 

Hotan. Bailes fc Ord I J2006h . 



van Straten 



d2003|); (4) 
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Table 2. PSR J0437-4715 timing model parameters' 1 



Parameter Name 


Parameter 


Tempo2 


Monte-Carlo 


Error 


and Units 


Value 


Error b 


Error b 


Ratio 




Fit and Data Set 


MJD range 


50191.0-53819.2 








Number of TOAs 


2847 








Rms timing residual (/is) 


0.199 








Measured Quantities 


Right ascension, a (J2000) 


04 h 37 m 15!8147635 


3 


29 


9.8 


Declination, 5 (J2000) 


-47°15 / 08 / /624170 


3 


34 


11 


Proper motion in a, fi a cos 5 (mas yr -1 ) 


121.453 


1 


10 


8.7 


Proper motion in 5, fi$ (mas yr -1 ) 


-71.457 


1 


12 


9.0 


Annual parallax, rr (mas) 


6.65 


7 


51 


7.9 


Dispersion measure, DM (cm -3 pc) . . . 


2.64476 


7 


d 




Pulse period, P (ms) 


5.757451924362137 


2 


99 


47 


Pulse period derivative, P (10~ 20 ) 


5.729370 


2 


9 


4.8 


Orbital period, Pb (days) 


5.74104646 c 


108 


200 


1.9 


Orbital period derivative, Pb (10~ 12 ) . . . 


3.73 


2 


6 


2.5 


Epoch of periastron passage, T (MJD) 


52009.852429 c 


582 


780 


1.3 


Projected semi-major axis, x (s) 


3.36669708 c 


11 


14 


1.4 


Longitude of periastron, u (°) 


1.2224 c 


365 


490 


1.3 


Orbital eccentricity, e (10~ 5 ) 


1.9180 


3 


7 


2.1 


Periastron advance, Co (° yr~ x ) 


0.01600 c 


430 


800 


1.8 


Companion mass, m 2 (M Q ) 


0.254 c 


14 


18 


1.3 


Longitude of ascension, Q (°) 


207.8 C 


23 


69 


3.0 


Orbital inclination, i (°) 


137.58 


6 


21 


3.7 



Set Quantities 
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Table 2 — Continued 



Parameter Name 


Parameter Tempo2 Monte-Carlo Error 


and Units 


Value Error b Error b Ratio 


Reference epoch for P, a 




and 5 determination (MJD) 


52005 


Reference epoch for DM 




determination (MJD) 


53211 



a These parameters are determined using Tempo2 which uses the In- 
ternational Celestial Reference System and Barycentric Coordinate Time. 
As a result this timing model must be modified before being used 
with an observing system tha t inpu ts Tempo format parameters. See 



Hobbs. Edwards fc Manchester I (120061 ) for more information. 

b Given uncertainties are la values in the last digits of the parameter values. 

c Because of large covariances, extra precision is given for selected parameters. 

d Dispersion measure was determined through alignment of simultaneous 
CPSR2 observations centred at 1341 MHz and 1405 MHz. The effect of red 
noise is therefore not applicable. 
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Table 3. Comparison of DE200 and DE405 results for PSR J0437-4715 



Parameter name 


DE200 result 


DE405 result 


Rms residual (ns) 


281 


199 


Relative \ 2 


2.01 


1.0 


Parallax, tt (mas) 


7.84(7) 


6.65 (7) 


Parallax distance, (pc) 


127.6(11) 


150.4(16) 


Previously published 7r (mas) 


7.19(14) 1 


6.3(2) 2 


Kinematic distance, Dk (pc) 


154.5 (10) 


156.0 (10) 


Dk corrected for Galactic effects (pc) 


155.5 (10) 


157.0 (10) 


Variation of Newton's gravitational 






constant, \G/G\ (KT 12 yr" 1 ) 


-21.2(22) b 


-5.0(26) b 


Total proper motion, \i tot (mas yr _1 ) 


140.852(1) 


140.915(1) 


Companion mass, m c (M Q ) 


0.263(14) 


0.254(14) 


Pulsar mass, m psr (M Q ) 


1.85(15) 


1.76(15) 


Periastron advance, Cj (° yr _1 ) 


0.020(4) 


0.016(4) 


GR prediction of Cj (° yr _1 ) 


0.0178(9) 


0.0172(9) 



a Numbers in parentheses represent the formal Tempo2 1 a uncertainty 
in the last digits quoted, unless otherwise stated. 

b Given are 2 a errors, i.e. 95% confidence levels. 



References. — (1) 

Hotan. Bailes fc Ord J2006h 



van Straten et al 



d2Q0l|); (2) 



